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Phase separation in the bosonic Hubbard model with ring exchange 
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We show that soft core bosons in two dimensions with a ring exchange term exhibit a tendency 
for phase separation. This observation suggests that the thermodynamic stability of normal bose 
liquid phases driven by ring exchange should be carefully examined. 
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Interest in ring exchange interactions in quantum 
many-body systems has a long history, both theoretically 
and experimentally 0. Recently, the ring exchange inter- 
action has been invoked in an effort to understand various 
aspects of high temperature superconductivity. While 
the Heisenberg model alone provides a rather accurate 
picture of magnetic excitations in the parent compounds 
of the cuprate superconductors 0, estimates of the mag- 
nitude of the ring exchange term are as high as one quar- 
ter of the exchange coupling[3j,y,|5( and it therefore has 
been of interest to understand how this term might mod- 
ify magnetic properties 0, IE IE El- Ring exchange in- 
teractions have also been suggested as a likely candidate 
to reconcile the properties of the underdoped pseudogap 
regime. The basic picture is that the ring exchange inter- 
action can give rise to a new normal "Bose metal" phase 
at zero temperature in which there are no broken symme- 
tries associated with superfluidity or charge density wave 
phases, and in which the compressibility is also finite |10|. 

With these motivations partly in mind, Sandvik et 
studied the phase diagram of the two-dimensional 
spin-1/2 XY model with spin exchange interaction on a 
square lattice, 



h = -jJ2b u -kJ2 Aiki 
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and (ij) denotes nearest neighbours and (ijkl) are sites at 
the corners of a plaquette. As is well known, for K = 
this model is exactly equivalent to the hard core bosonic 
Hubbard model, at half filling, with no interactions apart 
from the constraint on site occupations, and it has only a 
superfluid phase. Sandvik et al\l]$ studied the phase di- 
agram as a function of K using the Stochastic Series Ex- 
pansion algorithm They found that as K increases, 
the superfluid density, p s , decreases up to a critical value, 
K c , where a phase transition takes place, and p s goes to 



(tt,0) and (0,7r) channels of the plaquette-plaquette cor- 
relation function. This indicates the existence of a phase 
in which plaquettes with large and small values of ring 
exchange alternate with a striped pattern across the lat- 
tice. This phase is also an incompressible insulator. For 
yet larger K, charge density wave order is established 
in which the site occupations are alternatingly large and 
small. 

It is believed that when the hardcore constraint is 
relaxed, the striped plaquette phase might evolve into a 
normal compressible conducting "Bose metal" in which 
none of the order parameters mentioned above is non- 
zero. This suggestion leads us to study here the phase 
diagram of the soft core bosonic Hubbard model at half 
filling with ring exchange interaction, 

H = -i^afa, + a j ai ) + U^2,ni{m - 1) (3) 

<iJ> i 
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where the destruction and creation operators satisfy 
[ai,a|] = rii = a\a,i is the number operator at site 
i and U is the onsite interaction strength. For our quan- 
tum Monte Carlo (QMC) simulations we used the World 
Line algorithm with four-site decoupling|l7|. We verified 
our code for K = by comparing with existing results 
for hard and soft core bosons with and without near and 
next near neighbor interactions. For the K ^ case, we 
compared with the hard core results of 

Before discussing results for the full many-body sys- 
tem, it is interesting to study the behaviour of two 
bosons, since the formation of a bound state is closely re- 
lated to the issue of phase separation. Ring exchange, like 
an attractive potential, favors proximity of two bosons, 
since the action of such a term is nonzero only when 
the two bosons live on the same plaquette. In Fig. 1 we 
show the average separation (<& |r 2 |$o) °f tw0 bosons, 
normalized to the number of sites L 2 on an LxL lattice. 
As K increases, there is a crossover at K/t w 3 from a 
regime where the boson separation grows linearly with 
system size, so that the normalized separation is size in- 
dependent, to one in which the boson separation does not 
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FIG. 1: Average ground state separation ($o|r ,2 |$o) of two 
bosons, normalized to the system size, as a function of the 
magnitude of the ring exchange energy scale K. Here the hop- 
ping t = 1 and the soft-core repulsion U = 12. At small K, the 
normalized separation is independent of lattice size, indicat- 
ing that the two bosons are spread independently throughout 
the lattice. At larger K the bosons prefer smaller separation 
to optimize the ring exchange energy, indicating the formation 
of a bound state. 




While we do not show the associated data, plots of 
(<J>o |y 2 l^&o) for different soft-core repulsion U reveal that 
the average separation is insensitive to the value of U. 
This will obviously be true in the unbound regime at 
small K, since the density is so dilute. It is loss clear that 
this should be so in the bound regime at large K. How- 
ever, as can be seen from the data in Fig. 1, the radius of 
the bound state is several lattice spacings ((r 2 ) cx 0.1X 2 
whence r cx 0.3L), so here too the effect of U is expected 
to be relatively small. 

Fig. 1 suggests that there might be a tendency for ring 
exchange to cause the bosons to clump together, and, in 
an extreme scenario, to undergo phase separation. How- 
ever, at densities higher than the dilute two boson case, 
this effect is opposed by the repulsion U. The focus of 
this paper is to examine this competition and determine 
the phase diagram of the soft core case as a function of 
both U and K at half filling. 

The most straightforward indication of phase separa- 
tion comes from a real space image of the boson density 
during the course of a simulation. Fig. 2 shows the av- 
erage density distribution[13| for L = 16, U = 4 and 
K = 2.5. We see indications that the bosons undergo 
phase separation: At less than quarter filling (top panel) 
the bosons clump together into a compact region of high 
density. At densities above quarter filling, on a lattice 
with periodic boundary conditions, the number of occu- 
pied plaquettes is largest (and hence the ring exchange 
energy is most negative) for a configuration where the 
bosons stretch out in stripe across the lattice (bottom 
panel) . 



FIG. 2: Typical QMC results for the average density dis- 
tribution in the phase separated region. Here [7 = 4 and 
K = 2.5. Top panel: p=50/256; Bottom panel: p=128/256 
(half- filling). At half- filling, a stripe across our periodic 
boundary condition lattice maximizes the number of occupied 
plaquettes, and hence minimizes the ring exchange energy. 

acteristic of a large portion of the K—U phase diagram by 
examining the density-density correlation function and 
its associated structure factor, fixing U (or K) and scan- 
ning K (or U). As Fig. 2 illustrates, if the bosons phase 
separate, they may form a structure in which a set of 
contiguous sites of about half the system size will have 
appreciable boson occupation. The other half of the lat- 
tice is essentially empty. Therefore, if one examines the 
structure factor of the density-density correlation func- 
tion, 

S(fc x ,Ag = ^£C7(r)e- irk (4) 

r 

with 

C(r) = ^]>>(r>(r + r')), ( 5 ) 

r' 

one should observe a peak in S at small momentum, 
e.g. (2tt/£,0), (0,2tt/£) or (2n / L , 2n / L) depending on 
the precise orientation of the clump [lj|- By looking at 
the sum of the density structure factor at these three 
smallest momentum values we are sensitive to phase sep- 
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FIG. 3: The average structure factor, (S(2tt/L,0) + 
S(0,2n/L) + S{2n/L,2n/L))/3 versus K. We see a sharp 
increase in S at K = 2 for U/t = 4 and K w 4 for (7/i = 8. 
The simulations were done for /3 = 8, and p = 0.5. 



FIG. 4: The superfiuid density as a function of A' for U = 8. 
p s remains finite in the phase separated region, indicating 
that the bosons are delocalized across the clump of occupied 
sites. Inset: The hard core limit for which, instead, p a — > 0. 



roughly circular shape (Fig. 2a) or in some more elon- 
gated pattern (Fig. 2b). Fig. 3 clearly shows this be- 
haviour: For K < 2 at U = 4, S is very small at the 
relevant momenta. For K > 2, phase separation sets in. 
Data for 16 x 16 and 24 x 24 lattices are shown and their 
agreement indicates that this phase separation is not a 
finite lattice effect. The critical value of K grows roughly 
linearly with U . 

It is also interesting to understand the behaviour of 
the superfiuid density p s . One does not necessarily ex- 
pect p s to vanish when phase separation occurs. In fact, 
as is well known, p s ^ for the soft-core boson Hubbard 
model at all fillings, including commensurate density, if 
U/t is sufficiently small. Similarly, here it is possible that 
p s can survive in the dense region of the phase separated 
lattice[2(3, [2l|. Our simulations show that when phase 
separation first occurs, the populated region forms a band 
that spans the whole system. The bosons may then de- 
localize along that band, maintaining an (anisotropic) 
superfiuid density. We have found that when the bosons 
form such a band, the plaquette-plaquette structure fac- 
tor is also anisotropic and has long range correlations 
along the direction of the band. As K is increased fur- 
ther, the populated region of the lattice takes the form of 
an island. In such a case, the system may not be consid- 
ered a superfiuid in that one cannot establish superflow 
across the system. However, the bosons may still be de- 
localized over the extent of the island 16] . Fig. 4 shows 
p s versus K. 

By making scans like those of Fig. 3 at several values 
of U, we construct the phase diagram which we show in 
Fig. 5. Above the solid line, the system undergoes phase 
separation. 

How do these results connect with the previous studies 
in the hard-core limit? There we know that a phase tran- 
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FIG. 5: The phase diagram of the half- filled Bose Hubbard 
model in the U — K plane. Below the solid line, the system is 
superfiuid while above the line it phase separates. See text. 



plaquette phase [llj. We show in Fig. 6 the behaviour of 
the boson density p as a function of chemical potential p, 
in this hard-core limit. The jump in p across half-filling 
p = 0.5 shows that the plaquette ordered phase has a gap 
to the addition of bosons ( "charge excitations" ) , a result 
which is in agreement with Sandvik etal [13 ■ The slope 
of the p versus p is the compressibility k. Consequently, 
if this curve "bends backwards" , the system is thermody- 
namically unstable and undergoes phase separation[18|. 
While the data are not conclusive, we do see hints of an 
instability. For \p — 0.5 1 > 0.008, the slope is finite and 
corresponds to a normal superfiuid. For | p — 0.5 1 < 0.008, 
the slope is either very large, or perhaps negative, indica- 
tive of phase separation. 

So far we have addressed mostly the half filled case, 
and lower densities. It is of course of interest to examine 
higher fillings where the effect of U will be expected to 
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FIG. 6: The boson density, p, versus the chemical potential, 
p, for hard core bosons. The jump in p at p = i is associated 
with the nonzero gap in the phase with long range striped 
plaquette order. There are indications of a region of nega- 
tive compressibility k — dp/dp immediately adjacent to the 
gapped phase. 

for p = 1 and p — 1.5 and found in both cases that despite 
the higher densities, phase separation still sets in above a 
critical value K > K c of the ring exchange energy scale. 
This result is not so surprising since the on-site repulsion 
U and the ring exchange term both scale with density as 
p 2 . (Each of the four ring exchange creation/destruction 
operators picks up a factor of ps *Jn when acting on a 
site with occupation n.) 

In conclusion, we have shown that a sufficiently large 
ring exchange energy can lead to a thermodynamic insta- 
bility and phase separation. We determined the critical 
if as a function of the soft core repulsion U for a half- 



filled lattice and found roughly K c m U/2. We conclude 
that the soft core boson Hubbard model does not exhibit 
a normal Bose metal phase. The bosons are either super- 
fluid or undergo phase separation. This phase separation 
also takes place when the hopping parameter vanishes, 
t = 0, a limit examined for the quantum phase model in 
reference^^ but which did not find phase separation. 

Finally, let us comment on the implications of our work 
for the phase diagram of the spin- 1/2 quantum Heisen- 
berg model with a ring exchange term. The kinetic en- 
ergy term in the hard-core boson Hubbard model maps 
onto JJ2(SfSj + S?SJ) with exchange constant J = 2t. 
At the value U = 4t in our soft core model, double oc- 
cupancy is already very rare at half-filling, and hence we 
are almost in the hard-core limit. The value of the ring 
exchange energy scale required to drive phase separation 
for this U is K w 2t, or in other words, K w J. To repli- 
cate the near-neighbor coupling of the z components of 
spin in the Heisenberg model we must include a near- 
neighbor repulsion in the bose-Hubbard model, a term 
which clearly would suppress phase separation. Thus we 
expect ring exchange to have the potential to drive phase 
separation in the Heisenberg model only for K consider- 
ably greater than J. 

Note added: After the completion of this work 
a preprint by Melko, Sandvik and Scalapino 
(cond-mat/0311080) appeared where the phase dia- 
gram of the doped hardcore system with ring exchange 
was determined. 
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